Derivation of quantum work equalities using quantum Feynman-Kac formula 
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On the basis of a quantum mechanical analogue of the famous Feynman-Kac formula and the 
Kolmogorov picture, we present a novel method to derive nonequilibrium work equalities for isolated 
quantum systems, which include the Jarzynski equality and Bochkov-Kuzovlev equality. Compared 
with previous methods in the literature, our method shows higher similarity in form to that deriving 
the classical fluctuation relations, which would give important insight when exploring new quantum 
fluctuation relations. 
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Feynman-Kac (FK) formula originally found by Feynman in quantum mechanics [l[ and extended by Kac Q 
establishes an important connection between partial differential equations and classical stochastic processes. Briefly, 
assuming that in a continuous diffusion process the probability of a stochastic trajectory X started from a state x' at 
time t' is P[X\x' ,t']. The solution u(x',t') of the following partial differential equation 

d v u{x', if) = -C+(x', t')u(x', If) - g(x', t')u(x', t'), . , 

u(x',t' =t) =q(x') {L) 

has a concise path integral representation Q: 

.M^rf^feWH, (2) 

where C + is the Markovian generator of the diffusion. This is the famous FK formula in classical stochastic processes. 

Pioneered by Lebowitz and Sphon [H , this formula was also found very useful in studying fluctuation relations ■ 
In the past two decades, these important relations have greatly dee pen ed our understanding about the second law 
of thermodynamics and nonequilibrium physics of small systems (3, Il0l-f20| . Recently, finding quantum fluctuation 
relations is attracting intensive interests. Fruitful theoretical and experimental results [2lM30| have been obtained. 
To our best knowledge, however, there is no work explicitly using the FK formula. At first glance, the reason is 
very obvious, because the classical trajectory picture on which the FK is based is not available in quantum physics. 
Contrary to the intuition, in this Rapid Communication we use an isolated quantum system as an example to show 
that there indeed exists a quantum mechanical analogue of the classical FK formula and it is very useful to derive the 
quantum nonequilibrium work relations including the Jarzynski (2lL l23l . [24| and Bochkov-Kuzovlev equalities pfjj j . 

Kolmogorov picture and backward invariable. We start by introducing essential notations and a new picture that 
is a quantum-mechanical analogue of Kolmogorov's idea |31[ in classical stochastic theory. Although the picture 
is virtually equivalent to other well-known pictures, e.g. the Heisenberg picture, we will see later that it is very 
relevant with the time reversal concept. We assume that the closed quantum system is described by a time-dependent 
Hamiltonian H{t). The system's density operators p at two different times t and t' (< t) are connected by the 
time-evolution operator U(t'), i.e., 

P (t) = U(t)U\t')p(t')U(t')uHt). (3) 

Given an arbitrary observable F that does not depend explicitly on time, we define its Kolmogorov picture as 

F(t,t') = U(t')F*(t)U\if), (4) 

where the superscript H denotes the Heisenberg picture: F (t)=W (t)FU(t). On the basis of Eqs. (O and the 
expectation value (F}(t) at time t in the picture is 

Tr[Fp(t)]=Tr[F(t,lf)p(lf)}, (5) 
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and the equation of motion for F(t,t') with respect to t' is simply 

ihd t >F(t,t') = -[F(t,t'),H(t')) 



F{t,t' = t) = F. ^ 

We see that it is a terminal condition rather than initial condition problem. It is worth pointing out that Eq. (|6]) 
is very distinct from the motion equation of the same F(t,t') with respect to the forward time t if the Hamiltonian 
explicitly depends on time. 

Equation ([5]) has a trivial property: the time derivatives on both sides with respect to t' vanishes, or equivalently, 
Tr[F(t, t')p(t')] being a backward time invariable. The property is very analogous to that of the Chapman-Kolmogorov 
equation in the classical diffusion theory (32j . According to our previous experience which constructing a more general 
backward time invariable may lead into the classical fluctuation relations [9j, it would be very interesting to explore 
whether the same idea is still true here. Imitating Eq. (5) in Ref. Q, we find there is very analogous backward time 
invariable in the quantum case: 

Tr[Fp(t)]=Tr[F(t,t')p(t% (7) 

if the new operator F(t, t') satisfies 

ihd t ,F(t,t') = -[F(t,t'),H(t')]-F(t,t'){ind t ,p(t') + [p{t'),H{t')])-p-\t') 

+(\F(t,t'),A(t')]B(t') + F(t,t')[B(t'),A(t')])p- 1 (t') (8) 

and its terminal condition at t is assumed to be F, where the operators A{t'), B(t') and invertible density operator 
~p(t') are arbitrary. The proof of Eq. ([7]) is straightforward. The meaning of the last two terms on the right hand side 
will appear when we chooses ~p(t') to be the system's density operator p(t'), i.e. the term ifidt''p + [p,H] vanishing. 
The general Eq. (|SJ) seem uncommon in the quantum mechanics except a specific case: 

iW v F{t,t') = -[F{t,t'),H(t')]-F(t,t')0(t% 

F(t,t' = t) = F, W 
where 0(t') is an arbitrary operator. We may easily write down its solution given by 

F(t,t') = U{t')F H {t)T+e^ti dTU '^ ^ u ^U\t') (10) 
= U(t')F H (t)Q(t,t')U\t>), (11) 

where 7+ is the time-ordering operator. We simply name Eq. (| 10[) quantum FK formula because of its highly formal 
similarity to the classical FK formula ^j. However, we must remind the reader that the whole time-ordering term, 
which we specially denote it by a operator Q(t,t') for convenience, only indicates that the operator satisfies 

ihd t ,Q(t,t') = -Q(t,t')[U^{t')0(t')U(t')] (12) 

with a terminal condition Q(t,t'=t)=l. So far, we mainly concentrate on a formal development; the physical 
relevance of the quantum FK formula (|10p and the backward time invariable (O is not obvious. In the following we 
show that these results would lead into the quantum Jarzynski and Bochkov-Kuzovlev equalities if one chose specific 
p(t'), A(t'), and B(t'). 



Quantum Jarzynski equality. We assume that the closed quantum system is initially in thermal equilibrium with a 
density operator p cq (0)=e~^ H ^ / Z (0) , where the partition function ■2T(0)=Tr[e _/3ff ( * l ]=e _,3G ( ) , (3 is inverse temper- 
ature and G(0) is the initial free energy. At later times the system evolutes under the time-dependent Hamiltonian 
H(t). We choose ^4=5=0, p(t') to be the instant equilibrium state p q(i') =e_ ' 3 V^(0 with the instant partition 
function Z(t')=Tr[e- /3ff( *')]=e-' 3G '(*' ) . Equation ||SJ} then becomes 

ihd t ,F{t,t') = -[Ffat^Hit'y-ihFit^dt'PeqWpe^it'). (13) 

Obviously, the above equation follows the structure of Eq. ([9]) and especially 

0(t') = ih[d t ,e- fm (^e pH(t,) + pd t ,G{t')] (14) 

when we substitute the expression of p eq (t') into the "source" term of Eq. (Tl3|) . Intriguingly, in this case Eq. (TT2|) 
indeed has a very simple analytical solution 

Q(t,t') = [[/t(t)e-^ ff (')[/(t)]t/t(j' )e ^(t') [/ ( i ')]^[G'(t)-G( t ')]. (15) 
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Hence, on the basis of Eqs. (O, (fTO)) , and (fT5|) we obtain 

Tr[F Pcq (t)} = TrlC/^O^Wr+e^'^^^^^^^^'^^^tOPcq^')] (16) 
= Tv[F^{t)e~^ Hm ^e^ Ha ^U^(t')p COi {t')U{t')] e ^(t)-G( t ')]. (17) 

If F=l and i'=0 Eq. (|17|) is just the quantum Jarzynski equality on the inclusive work [24j : 

(e -^ H (t) e ^(0) )eq(o) = e -^AG(t) j (18) 

where AG(t)=G(t) — G(0), and we have used ( ) cq (0) to denote an average over the initial density operator p eq (0). 
Additionally, Eq. (fl7|) at t'=0 is also a specific case of the general functional relation given by Andrieux and Gaspard 
earlier; see Eq. (12) therein [26j . 

Bochkov-Kuzovlev equality. Here we consider a special realization of the time-dependent Hamiltonian [33| : a dy- 
namic perturbation H\(t) (t > 0) is applied on a system that is initially in thermal equilibrium with a time-independent 
.ffo, that is, the total Hamiltonian at later times is H p (t)=H +Hi(t). We have used the subscripts o and p to indicate 
"perturbed" and "original", respectively. Obviously, the system's initial density operator is p p (0)=p o =e~P H ° /Z Q , 
where the partition function Z =Tr[e-' 9H °]=e-' 3G °. Choosing H(t')=H , Atf^-H^t'), B(t')=p{t')=p in Eq. ©, 
we obtain the following equation 

iKd t ,F(t,t') = -[F(t,t'),H ] - QF(M'),ffi(*')] -mt'^mit')^- 1 

= -[F(M0,#p(O]-%O[Po,ffi(i')]Po 1 - (19) 

We have used a new symbol F to distinguish it from the previous F because they satisfy different equations. We see 
that the above equation also follows the structure of Eq. © , and especially 

0(f) = [e-' 3H °,iJ 1 (t')]e' 3 ^. (20) 

It would be interesting to check whether there is a simple analytical solution to Eq. (|12[) under this circumstance. We 
find that it indeed has: 

Q(t,t>) = lUl(t)e-P H °U p mul(t'y H °U p (t% (21) 
Using Eqs. (0, (fTU|) . and (f2Tj) we establish another equality given by 

Tr[F Po ] = Tr[t/ p (i0(^)p WT+e^^^'^^^^^^^'^^^^^^^OPo] (22) 
= Tr[(ir)H(t)e-^^)pWe^^)"(*')C/t( t ' )i0oC / p ( t ')]. (23) 



where (F)^ (t)=U^(t) FU p (t) . If F=l and t'=0 Eq. (|23|) is the quantum Bochkov-Kuzovlev equality on the exclusive 
work 

{e -mo)?(t) e PH o)o = 1 (24) 

that was proposed very recently in Ref. [1^], where ( ) indicates an average over the initial density operator p . 

On the other hand, a distinction between the original and perturbed systems is not absolute. In physics we may 
regard H p (t) as the original quantum system while H Q is perturbed system if the dynamic perturbation of the latter 
is thought to be —Hi(t). Exchanging the subscripts p and o and changing the sign before Hi into minus in Eq. (|19l) . 
we then obtain another equation given by 

ihd t ,F(t,t') =-[F(t,t'),H (t')}+F(t,t')[p p ,Hi(t , )]p- 1 . (25) 

Of course, the argument can be as well strictly proved by choosing H(t')=H p (t'), A(t')=Hi(t'), and B(t')=~p(t')=p p (t') 
in Eq. ©. We are interested in what new equalities like Eqs. (|2"2l and (f2"3")) will be yielded. Doing analogous derivations, 
we have the following results: 

Tr[Fp p (t)] = Tr^O^WT+e^ 1 ^ 1 # dTU ^P^- H ^ T ^ l ^ u -^U^(t')p p (t')] (26) 
= Tr[F I W(p p ) I (t)(p- 1 ) I (t')t/o t (<')Pp(i')t/o(i')], (27) 
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where F 1 ^) and p l p (t) arc the interaction pictures of the observable and density operator, respectively, e.g. 
F 1 (t)=Ul(t)FU {t). Particularly, if F=l and t'=0, Eq. flU} becomes 

((p P ) I We^)o = e^. (28) 

So far, we are not very clear what physics the equality reveals, though its first order approximation should be related 
to the famous fluctuation-dissipation theorems [33j . 



Time reversal. In the remaining part we want to give time reversal explanations of Eqs. (|13[) . (|19[) . and (|25p. which 
is essential to understand the physical meaning of the backward time t' and the origin of these three equations. We 
use Eq. (fT3")) as an illustration. Multiplying its both sides with p eq (t') and introducing a parameter s=t—t' (0<s<t), 
we rearrange the equation into 

ihd s [F(t,t-s)p cq {t-s)} = [F(t,t-s)p cq (t- S ),H(t-s)}. (29) 



Noting this is an initial condition rather than the terminal condition problem. Equation (|29|) seems very analogous 
to an evolution equation of a density operator, which is indeed true if we multiply both sides of the equation by the 
antiunitary time-reversal operator 6 and its conjugation to obtain 

ihd s p R ( S ) = {H K (s),p R (s)), (30) 

where the time-reversed density operator and time-reversed Hamiltonian are 

Tr[Fp cq {t)\ 

H R (s) = QH(t-s)Qt, (32) 

respectively, and the coefficient is for the normalization of p R (s). Equation pip immediately explains the backward 
time invariable in the Jarzynski equality case. Additionally, one may also prove that the same equation is equivalent 
to the key Lemma in Ref. [26j that was used to prove the functional relation; see the Appendix. Doing very similar 
calculations, we find that Eq. (|19[) is equivalent to the evolution equation of the time-reversed density operator Pr(s) 
with Hamiltonian H R (s), which are given by 

M») = 7 FT ^eF(t,t- s ) Po et, (33) 

H R (s) = GH p (t- s)9 t = Ho + QH^t- s)6 t , (34) 

respectively, and Eq. (|25j) is equivalent to another evolution equation, the time-reversed density operator and Hamil- 
tonian of which are 

Ms) = =-^eF(t,t-s)p p (t-s)Qi, (35) 

H R (s) = 6i/ o e t = H , (36) 

respectively. Comparing the above three time reversal explanations in the same dynamic perturbation problem, wc 
sec that, though the original and time- reversed Hamiltonian for the quantum Jarzynski equality is completely the 
same with that for the quantum Bochkov-Kuzovlev equality, the time-reversed density operators have very distinct 
initial conditions, i.e., ~p R (0)=p eq (t) and /5r(0)=/9 o . Intriguingly, for the equality (|28[) . the time-reversed Hamiltonian 
is overlapped with the original one, both of which are the unperturbed H a , while the initial conditions of the 
time-reversed density operator and the original one are different: they are p p (t) and p Q , respectively. 

Conclusions. In this work, we have used a quantum mechanical analogue of the classical FK formula to derive 
known quantum noncquilibrium work relations in isolated quantum systems. Compared with the previous methods 
in the literature, our method shows highly similar in form to that for the classical fluctuation relations which we 
developed earlier Q . We think that it is insightful because one may find the backward time invariable first and then 
give its physical interpretation rather than vise versa. Previous work has shown that direct defining nonequilibrium 
physical quantities was very nontrivial in quantum case [28j . Extending our method into more complicated quantum 
systems, e.g. the open quantum systems would be a challenge in our following researches. 



This work was supported in part by the National Science Foundation of China under Grant No. 11174025. 



I. APPENDIX 

Equation (|3ip could be further simplified, the time reversed density operator at later time s is connected with the 
initial condition by 

pr(s) = u R ( s )p R (o)ul( s ) = ( t)] u R (s)eF Pecl (t)e^u R (s), (37) 

where U R (s) is the time-evolution operator for the time-reversed Hamiltonian H R (s). Substituting the above equation 
and the solution of Eq. (fl~3"| 



F(t,t') = U(t')F H (t)e-^^e pHH ^U\t')e^- G ^ (38) 
into Eq. (|3f I) and doing a simple calculation we obtain 

U R (s) = eU{t~s)U\t)O f . (39) 
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